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Abstract 

We describe a semiclassical treatment of nuclear fusion reactions involving weakly bound nuclei. 
In this treatment, the complete fusion probabilities are approximated by products of two factors: 
a tunneling probability and the probability that the system is in its ground state at the strong 
absorption radius. We investigate the validity of the method in a schematic two-channel application, 
where the channels in the continuum are represented by a single resonant state. Comparisons with 
full coupled-channels calculations are performed. The agreement between semiclassical and quantal 
calculations is quite good, suggesting that the procedure may be extended to more sophisticated 
discretizations of the continuum. 

PACS numbers: 25.60.-t, 25.60.Pj, 24.10.-i 

Keywords: Unstable beams, Fusion reactions, Projectile breakup 



1 



I. INTRODUCTION 



The effects of channel coupling in fusion reactions induced by weakly bound projectiles 
have attracted great interest over the last decade ]| . Some theoretical studiespredict strong 
influence of the breakup channel over the complete fusion (CF) cross section (2, [J ^[ ^, fUll- 
One of the reasons for this interest is the growing amount of experimental data on these 
processes jsl S 10, [ill 12, 13, 3' which are important both for the study of astrophysical 
processes as well as for the production of superheavy nuclei. 

The appropriate theoretical tool to handle this problem is the coupled-channels method. 
However, its implementation becomes very complicated for the breakup channel, which is 
specially important in the case of weakly bound nuclei, since it involves an infinite number 
of states in the continuum. For practical purposes, it becomes necessary to approximate 
the continuum by a finite set of states, as in the Continuum Discretized Coupled-Channels 
method (CDCC) 1 151. This procedure has been extended to the case of fusion reactions 

n n □ 

in refs. p, la, 13]. Recently, a semiclassical alternative based on the classical trajectory 
approximation of Alder and Winther (AW) Q] has been proposed [l?]]. This approximation 
was used to calculate breakup cross sections and the results were compared with those of 
the CDCC method. The agreement between these calculations was very good. Since this 
semiclassical version of the CDCC method is much simpler, it may be a very useful tool to 
calculate cross sections for other channels in reactions with weakly bound nuclei. 

Although the AW method has been extensively used for several nuclear reaction pro- 

n 

cesses, only very recently it was applied to the estimate of the fusion cross section [18j. For 
this application it was considered a simplified two bound channels problem for which the 
fusion cross section obtained with the AW method was compared with results of a quantal 
coupled-channels calculation. In spite of the large simplification in the problem, the agree- 
ment between these two calculations at above-barrier energies was very good. However, the 
semiclassical method severely underestimated the fusion cross section at sub-barrier ener- 
gies. In the present work, we also consider two channels, but instead of being both bound, 
one of them represents a resonance in the continuum. We have assessed the validity of 
the semiclassical model in this case through comparisons with results of the corresponding 
full quantum mechanics coupled channels calculations. We have shown in particular that 
the good agreement with quantum mechanical calculations can be extended to sub-barrier 
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energies through the analytical continuation of the time. 

The paper is organized as follows: in section |H] we review the semiclassical treatment of 
nuclear reactions and give its extension to fusion reactions. In section 11111 we discuss the 
application of this treatment to the schematic two-channel model of ref. jl8j . and show that 
it severely underestimates the fusion cross section at sub-barrier energies. Then, in section 
IIV1 we develop a generalization of the semiclassical model which is accurate at sub-barrier 
energies and takes into account the width of the effective channel. The conclusions and 
discussion of further work are given in the last section. 

II. SEMICLASSICAL DESCRIPTION OF FUSION 

As this work is devoted to reactions induced by weakly bound projectiles, the variables 
employed to describe the collision are the projectile-target separation vector, r, and the 
relevant intrinsic degrees of freedom of the projectile, £. For simplicity, we neglect the 
internal structure of the target. The Hamiltonian is then given by 



where h (^) is the intrinsic Hamiltonian and V(r, £) represents the projectile-target interac- 
tion. The eigenvectors of /io(£) are given by the equation 



The AW method |l6( is implemented in two-steps. First, one employs classical mechanics for 
the time evolution of the variable r. The ensuing trajectory depends on the collision energy, 
E, and the angular momentum, hi. In its original version, an energy symmetrized Rutherford 
trajectory ri(t) was used. In our case, the trajectory is the solution of the classical equations 
of motion with the potential V(r) = (<p Q \ V(r,^) \cj) Q ) , where \4>o) is the ground state of the 
projectile. In this way, the coupling interaction becomes a time-dependent interaction in 
the £-space, Vi(^,t) = V(rj(t),£). The second step consists of treating the dynamics in 
the intrinsic space as a time-dependent quantum mechanics problem. Expanding the wave 
function in the basis of intrinsic eigenstates, 



h = h (0 + V(r,0 



(1) 




(2) 




(3) 



a 
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and inserting this expansion into the Schrodinger equation for i[>(£,t), one obtains the AW 
equations 

iha a (l,t) = Y. (MVi&Qlfo) e< e "- e ^ h ap(l,t). (4) 

P 

These equations are solved with the initial conditions a a (l,t — > — oo) = 5 a o, which means 
that before the collision (t — > — oo) the projectile was in its ground state. The final popu- 
lation of channel a in a collision with angular momentum I is P[ a ^ = \a a (l,t — > +00) | 2 and 
the angle-integrated cross section is 

° a = ^ £(2/ + l) P i a) - 

To extend this method to fusion reactions, we start with the quantum mechanical calcula- 
tion of the fusion cross section in a coupled channel problem. For simplicity, we assume that 
all channels are bound and have spin zero. The fusion cross section is a sum of contributions 
from each channel. Carrying out partial-wave expansions we get 
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(6) 



with 



Ah r 

P /»= ^ J dr \u al (k a ,r)\ 2 W?{r). (7) 

Above, u a i(k a ,r) represents the radial wave function for the /^-partial-wave in channel a 
and is the absolute value of the imaginary part of the optical potential associated to 
fusion in that channel. 

To use the AW method to evaluate the fusion cross section, we make the approximation 

P?( a )~P«>Ti a \E a ). (8) 

In eq.flBJ), is the probability that the system is in channel a at the point of closest ap- 
proach on the classical trajectory, and Ti a \E a ) is the probability that a particle with energy 
E a = E — e a and reduced mass \i = M P M T / (M P + M T ), where M P , M T are respectively 
the masses of the projectile and target, tunnels through the potential barrier in channel a. 

We now proceed to study the CF cross sections in reactions induced by weakly bound 
projectiles. For simplicity, we assume that the GS is the only bound state of the projectile 
and that the breakup process produces only two fragments, F 1 and F 2 . In this way, the labels 



a = and q ^ correspond respectively to the GS and the breakup states represented by 
two unbound fragments. Neglecting any sequential contribution, the CF can only arise from 
the elastic channel. In this way, the cross section <j C f can be obtained from eq.Q, dropping 
contributions from a^O. That is, 

^ = ^E( 2/ + 1 ) Pi ur "Tf\E\ (9) 
i 

where 

p5^p(0) = |ao(/;U| 2 (1Q) 

is usually called survival (to breakup) probability. 



III. COMPARISON WITH COUPLED-CHANNELS CALCULATIONS 

The accuracy of the semiclassical fusion cross section has recently been checked in a 
schematic two-channel calculation for the scattering of 6 He projectiles on a 238 U target at 
above-barrier energies ^|. The weakly bound 6 He nucleus dissociates into 4 He and two 
neutrons, with threshold energy B = 0.975 MeV. The elastic channel is strongly coupled to 
the breakup channel and the influence of this coupling on the fusion cross section is very 
important. In this calculation, the breakup channel was represented by a single effective 
state 3|. For simplicity, the effective channel was treated as a bound state but it was 
assumed to contribute only to incomplete fusion. Furthermore, the ingoing wave boundary 
condition is used in all these calculations. The CF cross section is therefore given by eq. 
(J2J). In ref. [3] the threshold energy was neglected and the same potential barrier was used 
for both channels, i.e. Wood-Saxon shapes for the real and imaginary potentials, with the 
parameters: V = -50 MeV, r 0r = 1.25 fm, a r = 0.65 fm, W = -50 MeV, r 0i = 1.0 fm 
and = 0.65 fm. This calculation adopted a form factor with the radial dependence of 
the electric dipole coupling and the strength was chosen arbitrarily, so that the coupling 
modified appreciably the cross section predicted by the one dimension barrier penetration 
model. The results (figure 1 of ref. are shown here in figure 1 (a). Above the Coulomb 
barrier, the semiclassically calculated cross section is in very good agreement with the one 
obtained by the coupled-channels method. However, the situation is much worse at energies 
below the Coulomb barrier. This is illustrated in figure 1 (b), where the results of (a) 
are plotted on a logarithmic scale. As it usually occurs, the total fusion cross section at 
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FIG. 1: Complete fusion cross section as calculated quantum mechanically (full line) compared 
with that of the Alder and Winther calculation (full dots) for the two channel case with B = 0. 
(See text for details). For comparison, the fusion cross section calculated by quantum mechanics 
in the no coupling case is also shown (dashed line). 

sub-barrier energies is enhanced in the CC calculation. The contribution from CF alone 
(solid line) is slightly larger than the fusion cross section in the no-coupling limit (dashed 
line). It is clear that the semiclassical calculation does not reproduce this effect, since its CF 
cross section is much lower in this energy range. This deficiency is due to the fact that in 
the semiclassical procedure described in section |H] the coupling is not considered inside the 
barrier region, because there is no classical trajectory connecting points inside this region. 
Therefore, the effective barrier lowering ocurring in the coupled-channels calculation cannot 
be reproduced. The semiclassical calculation shows also some instability for energies close 
to the barrier top (Vb — 20 MeV). This last behavior results from orbiting effects in the 
classical trajectories. 

IV. THE IMPROVED SEMICLASSICAL TWO-CHANNEL MODEL 

To improve the semiclassical model at sub-barrier energies, we resort to the analytical 
continuation method, which consists of introducing the imaginary part of the time variable 



6 




- oo 



wt=t R 



(B) t = t R 



(D) t = t R -iA 
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FIG. 2: Analytical continuation of the time variable. The upper panel shows the branches of the 
classical trajectory and the lower panel the evolution on the complex time plane. 



to obtain a classical trajectory in the sub-barrier region [20]. This procedure is illustrated 
in figure 121 where the time scale is chosen such that t — at the external turning point, 
r e . Along the incident branch of the trajectory (A), the time develops on the real axis as 
the system approaches the barrier. At r = r e , the trajectory splits into two parts: the 
reflected branch (B) and the classically forbidden transmission branch (C). On the former, 
which is not relevant to the fusion process, the time remains real. Along the branch (C), 
the real part of t remains equal to zero while its imaginary part develops on the negative 
part of the imaginary axis, until this trajectory reaches the exit point r iy at t = —iA. This 
trajectory is then continued into the internal classically allowed region (D), towards the 
strong absorption radius, Rp, where fusion occurs. Over this branch, the real part of t 
grows whereas its imaginary part keeps the value tj = —A. 

We used the analytical continuation of t in the calculation of the classical trajectories and 
generalized the AW equations (JH) accordingly. In figure 01 CF cross sections of different cal- 
culations at sub-barrier energies are compared. The open and the full circles are respectively 
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FIG. 3: Complete fusion cross section obtained with the calculations discussed in the text. 

results from semiclassical calculations with real and complex times. The solid line repre- 
sents results of coupled channel calculations. It is clear that the analytical continuation of 
the time variable improves substantially the accuracy of the semiclassical method. Now, 
their results agree very well with the ones of the quantum mechanical calculation, except 
for the last energy point, which is close to the barrier top. In this case, the cross section is 
strongly influenced by the orbiting effect in the classical trajectory. The system remains a 
very long time near the barrier radius, where the coupling is very strong. In this way, the 
elastic channel recovers the amplitude lost to the excited channel along the trajectory and 
the CF cross section becomes very large. This is a limitation of the semiclassical method. 
However, the situation tends do be much better in a realistic treatment of the continuum, 
with many channels and continuum- continuum coupling, owing to the irreversibility of the 
breakup process. The same occurs in the two-level model if the effective channel is treated 
as a doorway with finite width. This will be shown below. 

To further improve the semiclassical calculation, we take into account the reaction Q- 
value. We assume that the excitation energy of the effective channel is the breakup threshold, 
B. The next step is to simulate the irreversible nature of the breakup process, identified 
in the CDCC calculations of Diaz- Torres and Thompson [f| by an effective channel with 
complex energy. To this end we use the fact that an exponentially decaying state with mean 



8 




FIG. 4: Quantum calculation of the complete fusion cross section as a function of the cut-off radius 
of the potential —iT/2. For details see the text. 

life r = h/T, can be obtained through the inclusion of a constant imaginary potential equal 
to — iT/2 in the system Hamiltonian. This procedure, however, requires some care. The 
solution of the AW equations does not present difficulties since the population of the resonant 
state is vanishingly small as t — > — oo. The numerical solution of the coupled- channel 
equations, however, requires attention. To handle this situation one should switch-off the 
—iT/2 imaginary potential at some distance -R^ much larger than the range of the potentials, 
and then match the radial wave functions with their asymptotic forms. To illustrate this 
procedure, we show in figure 0] the CF cross section obtained through the solution of the 
coupled channel equations as a function of the cut-off radius. For this example, we considered 
the collision energy E c _ m , = 20 MeV ~ Vb and the typical width T = 2 MeV. The figure 
indicates that the result does not depend on the cut-off radius, provided that it is large 
enough (R^ > 20 fm). 
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FIG. 5: Complete fusion cross section as calculated quantum mechanically (full line) compared 
with that of the present work (full circles), including the analytic continuation in the classically 
forbidden regions, for the two channel case with B = 0.975 MeV and V = 2 MeV. 

Let us now consider the CF cross sections obtained with the above discussed procedures. 
The results of the improved semiclassical calculation (solid circles) are shown in figure in 
comparison with results of the CC method (solid line) and in the no-coupling limit. In order 
to exhibit the details above and below the barrier, the cross sections are represented on a 
linear (a) and on a logarithmic (b) scale. Firstly, one observes that the suppression of the 
CF cross section at above-barrier energies is less pronounced than in figure 1. The reason 
for this difference is that here we are taking into account the breakup threshold energy, 
B = 0.975 MeV, while the calculations of figure 1 were performed in the sudden limit 
(B = 0). Nevertheless, the suppression of the CF cross section remains quite. Comparing 
the semiclassical estimate for &cf with the CC values, we conclude that the improved 
semiclassical model leads to accurate results, above and below the Coulomb barrier. 

V. CONCLUSIONS 

We have illustrated how to extend the semiclassical method of Alder and Winther to the 
case of fusion reactions. Comparison of these calculations with full quantum ones for the 

10 



case of only two channels, the elastic and an effective one, shows a good agreement between 
them. Below the barrier we have used the analytic continuation of the trajectory in the 
classically forbidden region, and we have shown it essential to describe correctly the process. 

The fluctuations observed in the semiclassically calculated cross sections, in the case 
in which the effective channel has an infinite lifetime (zero width), may be associated to 
orbiting processes. However, when the resonance width is nonzero, the fluctuations vanish 
and the semiclassical calculations are in very good agreement with the full quantum ones. 
Thus it appears that fusion processes in weakly bound systems in which the elastic channel 
is coupled to the continuum are amenable to be treated by means of the extensions of the 
Alder and Winther's method presented here. 

The extension of this semiclassical method to less schematic treatments of the continuum 
is presently under way In it the continuum is represented by a discrete set of channels, 
as in |6j]; therefore the relative motion between the breakup fragments is more accurately 
described. In this way, detailed information on the positions and momenta of the projec- 
tile fragments at their creation point would be available. This would allow to study the 
collision of the fragments with the target, and therefore to estimate the incomplete fusion 
cross section as well as the sequential fusion contribution to the CF cross section. An esti- 
mate of this last process is not possible in the standard CDCC description of fusion [6j, so 
semiclassical calculations may furnish a way to solve this problem. It should be remarked 
that a semiclassical study al ong these lines has already been implemented for calculations 
of the breakup cross section [12| and the results were in very good agreement with the ones 
obtained with CDCC calculations. 

We acknowledge useful discussions with M.S. Hussein and P.R. Silveira Gomes. This work 
was supported in part by MCT/FINEP/CNPq(PRONEX) under contract no. 41.96.0886.00, 
PROSUL and FAPERJ (Brazil), and from PEDECIBA and CSIC (Uruguay). 
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